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Abstract: We apply a recent proposal for a distinguished ground state of a quantum 
field in a globally hyperbolic spacetime to the free massless scalar field in a causal diamond 
in two-dimensional Minkowski space. We investigate the two limits in which the Wightman 
function is evaluated (i) for pairs of points that lie in the centre of the diamond (i.e. far 
from the boundaries), and (ii) for pairs of points that arc close to the left or right corner. 
We find that in the centre, the Minkowski vacuum state is recovered, with a definite value of 
the infrared cutoff. Interestingly, the ground state is not the Rindler vacuum in the corner 
of the diamond, as might have been expected, but is instead the vacuum of a fiat space 
in the presence of a static mirror on that corner. We confirm these results by numerically 
evaluating the Wightman function of a massless scalar field on a causal set corresponding 
to the continuum causal diamond. 
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1 Introduction 

At the end of his classic textbook 'Aspects of Quantum Field Theory in Curved Space- 
Time,' [1] S.A. Fulling points out the central tension in the subject as traditionally studied: 

One of the striking things about the subject is the intertwining of the con- 
ceptual issues with the mathematical tools. Historically there has been a close 
association between relativity and differential geometry, while rigorous research 
in quantum theory has looked more toward functional analysis. Field quantiza- 
tion in a gravitational background brings these two alliances into head-to-head 
confrontation: A field is a function of a time and a space variable, 

x). 

Relativity and modern geometry persuade one with an almost religious intensity 
that these variables must be merged and submerged; the true domain of the 
field is a space-time manifold: 

ipix), X — ^ {x^} 



But quantum theory and its ally, analysis, are constantly pushing in the oppo- 
site direction. They want to think of the field as an element of some function 
space, depending on time as a parameter: 

Taking the relativistic view that the physical world has a space-time character indeed 
requires an approach to quantum field theory that is built on spacetime concepts, one that 
makes no fundamental appeal to "time as a parameter." The approach of algebraic quan- 
tum field theory takes this view seriously and attempts to push the foundations towards 
greater covariance, basing the theory on appropriate algebras of operators associated to 
open spacetime regions. The focus on the quasi-local algebras as the central constructs of 
the theory has encouraged the point of view that in general in curved spacetime there is 
no preferred quantum state and that a choice of state is akin to a choice of coordinates 
(see e.g. [2]). However, the algebraic approach has not been completely successful in 
constructing the expectation value of the stress-energy tensor or dealing with interacting 
fields. It remains an open question whether quantum field theory requires in addition the 
identification of a distinguished "ground state" (or class of them). Recently, a proposal 
has been made for exactly such a "ground state of a spacetime region" for a free quantum 
field theory [3, 4]. The proposal, implicit in work by S. Johnston on quantum field the- 
ory on a causal set [5], gives primacy to the 'true domain,' the spacetime manifold itself 
and its coordinate independent causal structure. The Wightman or two-point function, 
W{X, Y), of this "Sorkin- Johnston" (SJ) ground state is defined in an essentially covariant 
way starting from nothing more than A(X,Y), the Pauli- Jordan function or commutator 
function. Since the SJ ground state is taken to be Gaussian, its n-point functions are then 
completely determined in terms of W. 

In this paper, we study the SJ proposal in the case of a free massless scalar field in 
a causal diamond of two-dimensional Minkowski spacetime. We pay particular attention 
to the Wightman function obtained in the limit where the size of the diamond is large 
compared to the geodesic separation between the two arguments of W, and where the 
points lie either (i) far away from the left and right corners or (ii) close to one of the 
corners. 

In the former limit one might hope to recover the unique Poincare-invariant Minkowski 
vacuum state, if such a state existed. But in fact, no such vacuum exists, since W is 
Logarithmic and depends on an arbitrary length-parameter or "infrared cutoff" , as is well 
known. In our finite diamond, we find that W has the expected form, but with a definite 
value of the length-parameter determined by the diamond's area. 

In the latter limit, one might expect to see the Rindler vacuum state, since the geometry 
of the corner approaches that of the familiar Rindler wedge as the diamond size tends to 
infinity. However, we find that this is not the case. Instead, the SJ state close to the corner 
is the vacuum state of Minkowski spacetime with a static mirror on the corner. 

Further, we use the original causal set QFT formulation [5] to construct the ground 
state on a causal set that approximates the continuum causal diamond. We compare the 
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results with the foregoing continuum analysis in the subregions of the causal set corre- 
sponding to the two limits (i) and (ii). In both cases, the Wightman function on the causal 
set is in good agreement with the continuum Wightman function. 

2 The SJ vacuum 

We state the SJ proposal for the "ground state" of a free scalar field in a d-dimensional 
globally hyperbolic spacetime {M.,gfj_i,) [3, 5]. The starting point is the Pauli- Jordan func- 
tion A{X,X'), where we use capital letters X,X' to denote spacetime points. It is defined 

by 

A{X, X') := Gr{X, X') - Gr{X', X) (2.1) 

where Gr is the retarded propagator, i.e. the Green function of the Klein-Gordon operator^ 
□— satisfying the retarded boundary condition, Gr{X^ X') = unless X' ~< X, meaning 
that X' is to the causal past of X. Notice that A is real and antisymmetric. 

The rules of canonical quantisation relate the Pauli-Jordan function to the commutator 
of field operators by: 

iA{X,X') = [4>{X)J{X')] , (2.2) 
while the Wightman two-point function of a state |0) is 

Wo{X, X') = (O|0(X),^(X')|O) . (2.3) 

The Wightman function of the SJ state is defined by the three conditions [3]: 

1. commutator: iA{X, X') = W{X, X') - W*{X, X') 

2. postttvtty: dV dV nX)W{X, X')f{X') > 

3. orthogonal supports: fj^ dV W{X, X') W{X', X")* = 0, 

where J dV = J d-^X -g{X). On this view the Wightman function is a positive bilinear 
form on the vector space of complex functions. The first two conditions must be satisfied 
by the two-point function of any state. It is the third and last requirement that acts as 
the "ground state condition." 

For the conditions to specify a state fully, their solution must be unique. Suppose 
Wi and W2 both satisfy the conditions. Let us express the ground state condition as 
WW* = 0^. We have WiW^ = W^Wi = W2W^ = W^W2 = and Wi - W^ = W2 - W^. 
This implies 

{Wi - Wlf = {W2 - W2f 
(Wi + W^f = (W2 + w^f . 

^The signature of the spacetime metric is ( — h + ■ • ■ +) and □ = g'"'V^Vj/. 

^Strictly speaking, to multiply two quadratic forms together requires a metric, here given by a delta 
function. In these formulas the star denotes numerical complex-conjugate, not hermitian adjoint. 



-3- 



Both Wi + Wi and W2 + W2 are positive bilinear forms, and if there are unique positive 
square roots of such forms then this implies 

Wi + = W2 + W2 

=^ Wl = W2. 

Formally, we can interpret the ground state condition as the requirement that W be 
the positive part of iA, thought of as an operator on the Hilbert space of square integrable 
functions L'^{A4,dV) [3]. This allows us to describe a more or less direct construction of 
W from the Pauli- Jordan function [4-6]. The distribution iA{X,X') defines the kernel 
of a (Hermitian) integral operator (which we may call the Pauli-Jordan operator lA) on 
{M. , dV) . The inner product on this space is 

{f,g) := [ dVfiXygiX), (2.4) 
JR 

where dV = \/ —g{X)d'^X is the invariant volume-element on M.. The action of iA : / 1— >• 
iA/ is then given by 

{iAf){X):= f dV'iAiX,X')f{X'). (2.5) 
JM 

When this operator is self-adjoint, it admits a unique spectral decomposition. Whether 
or not this (or some suitable generalization of it) is the case depends on the functional 
form of the kernel iA(X, X') and on the geometry of A^. For the massless scalar field on a 
bounded region of Minkowski space, such as the finite causal diamond which is the subject 
of this paper, iA is indeed a self-adjoint operator, since the kernel iA(X, X') is Hermitian 
and bounded (see (4.1) below). In the following we assume that we can expand iA in terms 
of its eigenfunctions. 

Noting that the kernel iA is skew-symmetric, we find that the eigenfunctions in the 
image of iA come in complex conjugate pairs and T~ with real eigenvalues ztA^: 

(iAr±)(x) = ±A,r±(x), (2.6) 

where (X) = [T^{X)]* and Xq > 0. Now by the definition of iA, these eigenfunctions 
must be solutions to the homogeneous Klein-Gordon equation. If they are L^-normalised 
so that ||T+|p := (T+,T+) = l,"' then the spectral decomposition of the Pauli-Jordan 
function implies that the kernel can be written as 

iA{X,X') = Y,\T+{X)T^{X'r -Y,\T^iX)Tq{Xr- (2.7) 

We now construct the SJ two-point function Wsj{X,X') by restricting (2.7) to its positive 
part: 

Wsj{X,X') :=J2\T+{X)T+{Xr = r,(X)r,(X')^ (2.8) 
^ When the spacetime region is non-compact this should be replaced by a delta-function normalisation 

(r+,T+) = %-Q'). 
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where Tq{X) ■.= T+{X)^,. 

Let us compare this with the famihar case of the vacuum state of a free scalar field 
when the spacetime admits a timelike Killing vector k = dt- The free field can be expanded 
in a basis of modes, {u^, u^}, which are positive and negative frequency with respect to k: 

Uk{X) = Uk{t, X) = iVfce-*"^*?7fe(X) , (2.9) 

where Nj. is a normalisation- factor, uj. > 0, and X are the spatial coordinates of X. When 
the emission and absorption operators that serve as the expansion coefficients are given 
their customary normalisation, then the Pauli-Jordan function can be expressed as a mode 
sum 

zA{X, X') = [uk{X)ul{X') - ul{X)uk{X')\ . (2.10) 

k 

The vacuum state corresponding to the choice of positive frequency modes Uk is the ground 
state of the Hamiltonian associated to k, and its Wightman function is 



W{X,X') = Y,Uk{X)ul{X'). (2.11) 



Comparing (2.11) with (2.8) we see the eigenfunctions Tq := yJ~\qT^ chosen by the SJ 
condition are the analogues of the positive frequency modes in the static spacetime, but 
the SJ construction does not require the existence of a timelike isometry. 

When there is a timelike Killing vector, we can show that the SJ state (formally ex- 
tended to the case of a spacetime of infinite volume) is the ground state of the Hamiltonian 
associated with that Killing vector by showing that the Wightman function defined by 
(2.11) satisfies the SJ conditions. Clearly, (2.11) satisfies conditions (i) and (ii), and it only 
remains to check the ground state condition (iii). We have 

WW*{X,X") = [ dV'W{X,X')W*{X',X") 

Jm 

= j;^nfc(X)uKX") / dtd'^-'X'^^^)uUX')unX') 

k I '' 

= 5^5^nfc(X)nKX") j d"-iX'y=^C7,*(X')^r(X') |_ 



oo 
oo 



= EE^'^(^)^'(^") / ^^"~'X'v/^^C/,*(X')^;(X') / dte-^(-^+-')* 

since the sum over modes does not include the zero-mode ojk = 0- Thus, if the SJ conditions 
specify a unique state, that state is the appropriate ground state when there is a globally 
timelike Killing vector. 

This might be a good place to comment on the question of whether the SJ vacuum 
obeys the so-called Hadamard condition on its short-distance behavior. In static spacetimes 
the Hamiltonian vacuum obeys this condition, so the SJ vacuum does as well, as we have 
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just seen. On the other hand, Fewster and Verch [7] have recently provided examples of 
regions where the condition does not hold. In our opinion, the significance of the Hadamard 
condition will not be known until we understand better the nature of quantum gravity 
and its semiclassical approximations. Outside of that context Hadamard behavior seems 
irrelevant "operationally", since it corresponds in the Wightman function to the absence 
of a term that could only be noticed at extremely high energies. We hope to return to this 
question in another place, and to show in particular how, by "smoothing the boundary", 
one can tweak the ultraviolet behavior of the SJ state so that it becomes Hadamard. 



3 The massless scalar field in two-dimensional fiat spacetimes 

As background for our investigation of the massless scalar field in a two-dimensional causal 
diamond, we review the massless scalar field in two-dimensional Minkowski and Rindler 
spacetimes. The metric on Minkowski spacetime in Cartesian coordinates (t, x) is given by 

ds\j = -dt^ + dx^ . (3.1) 

Since the spacetime is globally hyperbolic and static with timelike Killing vector km = dt, 
we can separate the solutions to the Klein- Gordon equation into positive and negative 
frequency with respect to km- The field equation is 

UM(t> = -d'f^ + dlcp = ^ (3.2) 

and the normalised positive frequency modes can be taken as 

uf{t,x) = -^e-'^^'+'^'', (3.3) 

where uj^ = \k\. 

If we try to define a vacuum state, |0j\/) in the usual way as the state annihilated 
by the operator coefficients of the positive frequency modes in the expansion of the field 
operator (p{t,x), then it is well-known that we encounter an infrared divergence [8-10]: 

WM{t,x;t',x') := {OM\kt,x)^it',x')\OM) 

= Jl /"°° ^ -i\k\{t-t')+ik{x-x') (^-^^ 
4vr \k\ 

which is logarithmically divergent at A; = 0. 

Following [9], we can remove the divergence by introducing an infrared momentum 
cutoff A 



1 f°° dk 

— lim 

47r £-^0+ 



dA; [ ik{At+Ax-ie) I -ik{At-Ax-ie) 



lim [In [i(At + Ax - ie)fj] + In [i(At -Ax- ie)fi]] + 0(\A) 

47r e-s>0+ 

~lnn\d\ - Ugn{At)e{At^ - Ax^) + 0{XA), 
2tt 4 



(3.5) 
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where ^ = Xe"', 7 is the Euler-Mascheroni constant, and At = t — t' , Ax = x — x' , and 
d = \/—At'^ + Ax'^. The logarithm here is given a branch cut on the negative real axis 
and In denotes its principal value. The quantity A here stands collectively for At and Ax, 
such that small AA implies that both coordinate distances At and Ax are small compared 
to With non-zero A, the theory has a preferred frame. However, if we drop the 0{X) 
term in (3.5), we obtain a one-parameter family of two-point functions that depend on A 
but are fully frame-independent: 

WMxit, x; t', x') := -^lnfi\d\ - ]sgn{At)9{Af - Ax^) (3.6) 

Unfortunately, (3.6) cannot itself serve as a physical Wightman function, because it fails 
to be positive as a quadratic form (this being condition (ii) above). Nevertheless we will 
see that the form (3.6) will emerge in a natural manner as a certain limit of the two-point 
function we will derive for the diamond.^ 

It is worth noting that the theory whose fundamental field is the gradient of (j) rather 
than (p itself is not infrared divergent, and in fact the vacuum expectation value 

(OAf|V^<^(t,x)<^(t',x')|OA/) = 2n{AT-Ax^) ^^'^'^ 

already converges, except for the singularity on the lightcone. 

The Rindler metric [11, 12] arises from the Minkowski metric via the coordinate trans- 
formations t = a~ ^ e"'^ sinh arj and x = o~^e"^ cosha??, where a > is a constant with 
dimensions of inverse length and —00 < ^,ry < 00. The coordinates ^ and rj only cover a 
submanifold of the full Minkowski space, namely the right Rindler wedge, x > |t| ; but this 
submanifold is conformal to all of Minkowski space as one sees from the form of the line 
element in ^-r] coordinates: 

dsl = (-(i7?2 + ^^2) . (3.8) 

Lines of constant ^ are hyperbolae that correspond to the trajectories of observers accel- 
erating eternally at a constant acceleration ae~°'^ (figure 1), and are integral curves of the 
Killing vector = 9^. 

Since Rindler spacetime is globally hyperbolic and static in its own right, the canonical 
quantisation of the scalar field can be carried out in a completely self-contained manner [13]. 
Thanks to the conformal invariance of the massless theory in two dimensions, the field 
equation D^j^ = in Rindler coordinates (3.8) is just the usual wave equation, whose 
normalised positive frequency solutions, the FuUing-Rindler modes, are given by the plane 
waves 

= -^=e-^^v^+ivi ^ (3.9) 
y^47ra;p 

where Wp = Trying to define a Fulling- Rindler vacuum state |0/j) as the state annihi- 
lated by the operator coefficients of the positive frequency Rindler modes in the expansion 



^Perhaps (3.6) could also be understood as defining an "approximate state" valid when C\t and Ax are 
small compared to the IR scale A~^. 
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Figure 1. The right and left Rindler wedges of two-dimensional Minkowski spacetime. 



of (j) once again results in an infrared divergent integral for the two-point function W . Pro- 
ceeding as before, one can introduce a cutoff at small "momentum" A and discard a term 
of order A to obtain a one-parameter family of two-point functions, which are the same 
functions of Rindler coordinates as the Minkowski two-point function (3.6) is of Cartesian 
coordinates: 

W^i?,A(r/,^;r?',e') := -^In/i^lAr/^ - A^^, _ lsgn(A?7)e(Ar?2 - A^^) (3.10) 

where // = . As before, this two-point function is symmetric (boost-invariant) but fails 
to be positive and depends explicitly on the cutoff A. At best it can have an approximate 
validity when the coordinate-differences Ar/ and A^ are small compared to A~^. 

It is well known that the Minkowski and Fulling- Rindler "ground states" corresponding 
to (3.6) and (3.10) are not equal, since the Rindler mode functions are linear combi- 
nations of both positive and negative frequency Minkowski mode functions and u^* . 
This phenomenon is by now well understood as an instance of the Fulling-Davies-Unruh 
effect [13-15]: if the Rindler wedge is understood as a subregion of Minkowski space, and 
the field is in the usual Poincare-invariant vacuum state (say in 3 + 1 dimensions, where 
the latter is well defined), observers that are confined to the wedge and that accelerate 
eternally at a uniform rate will feel themselves immersed in a thermal bath of particles. 

The preceding calculations suffer from the appearance of infinite integrals and the 
consequent need for infrared cutoffs. Nevertheless, we will see that the two-point functions 
we have obtained in this section can be related to the case we study next, that of a finite 
two-dimensional causal diamond, where the the SJ construction and the integrals it gives 
rise to are completely well-defined. In this connection we comment also that inasmuch 
as both the Minkowski and Rindler spacetimes possess globally timelike Killing vectors, 
the formal demonstration of section 2 would apply to show that the SJ vacua of these 
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Figure 2. The causal diamond, u,v G {—L,L). The shaded portions of the diagram represent the 
centre («) and corner (ii) regions of interest in section 4.1. 

two spacetimes would coincide with the ground-states discussed in this section - were such 
states actuahy to exist. (See also [4].) 

4 The massless SJ two-point function in the flat causal diamond 

A causal diamond (or Alexandrov open set) is the intersection of the chronological future 
of a point p with the chronological past of a point q>~ p. Because such a causal diamond is 
a globally hyperbolic manifold in its own right, the scalar field possesses therein a unique 
Pauli- Jordan function A. In this section we will follow the SJ procedure to derive from A a 
two-point function W for the massless scalar field in a causal diamond in two-dimensional 
Minkowski space. We will then analyse the limit of W for points (i) in the centre of the 
diamond and (ii) in the corner of the diamond. 

Lightcone coordinates u = {t + x)/\/2 and v = {t — x)/\/2 are most convenient in this 
section. A causal diamond centred at the origin u = v = that corresponds to the region 
u,v £ {—L,L) is shown in figure 2 and will be denoted by Cl. Its spacetime volume is 
V = AL? . The commutator function il^{X,X') can be calculated in lightcone coordinates 
from the retarded and advanced propagators [16] : 



Given that Cl is a bounded region of spacetime, the associated integral operator iA. defined 
in (2.5) is of Hilbert- Schmidt type, i.e. its Hilbert-Schmidt norm is finite: 



iA(M, v; u , v') = — - [0(ii — u ) + 9{v — v') — 1] . 



(4.1) 



•L 

du dv du / dv^\iA{u^v]u\v^)\ = 2L^ < oo. 



(4.2) 



Since iA(X, X') = [iA{X' , X)]* , iA is also self-adjoint and so the spectral theorem applies. 
The positive eigenfunctions that satisfy iAT^ = +AfcT^, are given by the two sets [6] 



iku 



fk{u,v) := e 
gk{u,v) := e-'^"" 



+ e 



-ikv 



-ikv 



with k 



mr 



n 



1,2,... 



2 cos(fcL), 



with G /C, 



(4.3) 
(4.4) 



where /C = {A; G R | tan(A;L) = 2kL and A; > 0} and their eigenvalues are = L/k.^ Their 
L^-norms are ||/fc|P = SL"^ and = 8L^ — 16L^cos'^(A;L). It is clear that the eigenfunc- 

tions Tj^ with negative eigenvalues are given by the complex conjugates = [T^]*. To 
verify that the functions (4.3-4.4) and their complex conjugates are indeed all the eigen- 
functions with non-zero eigenvalue, we can use the fact [6, 17] that the sum over the squared 
eigenvalues of a Hilbert-Schmidt operator is equal to its Hilbert-Schmidt norm. A short 
calculation shows that 



^ 2L4 
^ (vrn)2 

n=l ^ ' 



2L^ 



2L^ lOL"^ 
~6~ ^ 6 



2L^ 



(4.5) 



coincides with (4.2), as required (the analytic evaluation of the second sum in (4.5) can be 
found in [6, 18]). The SJ prescription defines the two-point function Ws,j,l{u,v;u\v') as 
the positive spectral projection of iA: 



Wsj,l{u,v]u\v') 



n 



1 Trn ll/fel 



Jkiu,v)fkiu,v') + 



L 1 



9kiu,v)gl{u,v'). (4.6) 



We denote the two sums in (4.6) by 5i and 52, respectively. 
The first sum 



1 ^ -I 

1 V — > 1 r iunn 

Si = — > - e L 

n=l 



ivmr 

e L 



e L 



e L 



(4.7) 



can be evaluated in closed form. We recognise four Newton-Mercator series which converge 
to the principal branch of the complex logarithm: 



1 - e L 



In 



ITT (u — ) 

1 - e i — 



+ln 



1 - e E 



+ In 



i-K{v~u ) 

1 - e E 



(4.8) 



The second sum is 



^—iku _j_ ^—ikv 



2 cos{kL) 



^iku _j_ ^ikv 



2 cos{kL) 



keic 



kL [8 - 16cos2(A;L)] 



(4.9) 



We have no closed form expression for this sum but as n — t- oo, the roots of the transcen- 
dental equation tan(3;) = 2x rapidly approach Xn = ^^"2^^^ with n G N (see figure 3). 



^ There is a sign error in the commutator function in [6] which results in the fk and gt given here being 
the complex conjugates of those there. 
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Figure 3. Plot of tanx and 2x along with vertical lines at ^"^ ^ tt for n S Z. The values of the 
summation variable are the positions of the intersections tana; = 2a; > 0. 



Therefore if we approximate the sum by replacing k £ IC with Xn/L, we can expect the 
main error to come from a few modes of long wavelength. Consequently we can expect the 
resulting error term to be a slowly varying and small correction to the approximated sum. 
Let /Co := { ^^^"^ | n = 1, 2, 3 . . .} and define the error e{u, v; u', v') by 



5,= ^ 
fce/Co 



ZL-TT ^ 



2cos{kL) 



kL [8- 16cos2(A;L)] 



+ e{u, v; u' , v' 



47r ^ 2n - 1 

n=l 



iu(2n — l)7r 



iu(2n — l)7r 
2L -)- e 21 



iu\2n—l)TT iv\2n—l)TT 

e 2i + e ^ 



+ e{u,v;u',v') . 

(4.10) 

The sums above converge to logarithmic terms when L < oo and are most conveniently 

expressed in the form: 



S2 



1 



tanh ^ 



Z7r(Ti — u') 

e 2X 



+ tanh"^ 



+ tanh~^ 

i7r{u — v') 



i-Kiy — v') 

e 2X 



+ tanh"^ 



i-R^v—u) 

e 2Z 



(4.11) 



+ e{u, v\ u' , v' 



The two-point function of the SJ ground state in the causal diamond is given by the 
sum, Wsj,L = S\ + 82- Using tanh~^(x) = |ln(l + x) — |ln(l — x) and the fact that 
ln(l — e^) — ln(l + e^) = ln(l — eS), we can combine the two sums to yield 



-In 



1 — e 2X 



In 



1 — e 2X 



+ln 



iir{u—v') 

1 — e 2Z 



+ ln 



iir{v—u') 

1 — e 2Z 



+ e{u, v; u', v') 



(4.12) 



Wbox,L + e{u,v;u',v'), 



where Wbox,L is the exact continuum two-point function of the ground state of a massless 
scalar field in a box with reflecting boundaries at a; = ±^/2L. Wc shall now investigate the 
form of the SJ ground state in the limits (i) and (ii) mentioned above. 
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4.1 The SJ state in the centre and corner 

The limits we are concerned with require that two spacetime points be separated by a 
smah geodesic distance compared to the diamond scale L and that they be confined to (i) 
the centre of the diamond or (ii) a region near the left/right corner (we choose the left 
corner without loss of generality). In these particular limits, we want the arguments in the 
exponentials of the sums above to have a small magnitude so that we can Taylor expand the 
functions and obtain more illuminating forms of the two-point function. Keeping in mind 
that in u, ^-coordinates, the centre of the diamond lies at (0, 0), the limit that corresponds 
to a pair of points in the centre {i) can be defined as 

\u — u\<^L , \v — v'\<^L, \u — v'\<^L, \v — u'\ <^ L. (4.13) 

The limit that corresponds to the corner can be obtained in a similar manner by first 
translating the coordinate system such that the left corner of the diamond lies at the 
origin (0, 0) of the new coordinates. This corresponds to the (passive) coordinate transfor- 
mation X — )• X — V^L, or M — )• n — L and v — )• v + L. By "in the corner" we then mean the 
limit (a) where we first perform this translation and then apply the restriction (4.13) to 
the translated coordinates. 



The inequalities (4.13) constrain the pairs of spacetime points to be separated by 



a small geodesic distance \d\ = {2\u 



u V 



V 



<^ L and furthermore imply that 



\x'\ <^ L. This means that the points are confined to a narrow vertical strip centred 
on (i) the centre of the diamond or {ii) the left corner of the diamond, as illustrated in 
figure 2. Let the width of the strip he D <^ L. 



4.1.1 The centre 

We first analyse the sum in the centre by expanding to lowest order in 6/L, where S 



collectively denotes the coordinate differences u - 
ln(x) + 0{x) we identify the leading term in Si as 



-u'. Using ln(l — 



where Ci = 
tanh~^(e^l 



1 

8^ 



-lulu 



U V 



v'\ + ln|n 



V V 



u'\+Ci 



ITT 



+ 



(4.14) 



sgn(u — u') + sgn(t' — v') — sgn(u - 
= -glnf + 0{x) to obtain 



-sgT].{v — u'). Similarly, in 82-, we expand 



Si 



1 

8^ 



-ln|M 



u \v 



V 



ln\u 



V \\v 



41n h C2 — 

4L 2 



+ e + 0{-], (4.15) 



where C2 = sgn(u — u') + sgn{v — v') + sgn(u — v') + sgn{v — u'). 

Now we deal with the correction, e. Over a small region, e should not vary much. To 
investigate this, we now further restrict the arguments of W to lie in a small square centred 
on the origin, of linear dimension D, the width of the strip. After an analysis and numerical 
investigation given in the appendix, we find that e is indeed approximately constant over 
the small diamond and tends to a value ecentre ^ —0.063 as L tends to infinity. 
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The terms with arguments u — v' and v — u' cancel in Si + S2 so we obtain: 

Wcentre{u,v;u',v') = - -^ln| AnA?;| - ^sgn( An+ Ai;)6'( AtiAu) - ^In^ +ecentre +C ^-^^ , 

(4.16) 

as L gets large. Recall now that |AtiAf;| = It is then evident that (4.16) matches the 

"cut off" Minkowski two-point function Wm,x (3.6) with a particular value of the cutoff A 
given by 

A = ^ exp(-7 - 27reccntrc)^~^ ~ 0.46 x L^^ , (4.17) 

where 7 is the Euler-Mascheroni constant. As one would expect, A~^ ~ L for large L, 
leading to a logarithmic factor of the form ln|L^/AttAf |. 

Whilst strictly speaking we cannot take the L — t- 00 limit of such an expression, it 
seems fair to say that the SJ state takes on the character of a Minkowski vacuum in the 
centre of a large diamond. Notice finally that in the limit, the imaginary part of the 
two-point function does satisfy the requirement 9(VF) = A/2, as had to be the case. 

4.1.2 The corner 

For spacetime points close to the edges of the diamond, the boundaries of the diamond 
appear like the causal horizons of a Rindler wedge (see figure 1). To evaluate the SJ two- 
point function in this limit, we perform the translation u — )• n — L and v ^ v + L described 
above, which shifts the corner of the diamond to the position (0, 0). We then Taylor expand 
using (4.13) in the new coordinates. 

The sums Si and S2 can be evaluated as before, but the translation introduces it signs 
multiplying integer multiples of in. A brief inspection shows that the first sum (4.7) is 
unaltered by the translation, since only factors of exp(27ri) arise, while the second sum (4.9) 
picks up factors of exp(7ri) that induce sign changes in the terms involving u — v' and v — u'. 
The second sum now evaluates to 
1 



-ln|M — u'\\v — v'\ + ln\u — v'\\v — u'\ + C2 — 



+ e + 0{^], (4.18) 



where C2 = sgn(n — u') + sgn{v — v') — sgn(u — v') — sgn(u — u') = Ci. The correction 
term e can again be analysed numerically — see appendix A — and the result is that it 
varies very little over the small corner region for fixed L and tends to zero as L — )■ cxd. A 
consequence of the sign changes is that the constant terms that depend on L cancel in S2, 
whence there is no longer any obstruction sending L — t- 00. Taking this limit, we obtain 
the two-point function 



lim WcoTneT{u,v]u\v') = - -"-In 
L-!>oo 47r 



{u — u'){v — v') 



{u — v'){v — u') 



(4.19) 



^sgn(AM + I\v)[e{AuAv) - e{{v' - u){u - v))], 



which can be recognised as the two-point function of the scalar field in Minkowski space 
with a mirror at rest at the corner x = {x = —\/2L in the original coordinates) [19, 20]: 

Wcorner{t,X;t',x') = WM,x{t, X] t' , x') - WM,x{t, X] t' , -x') . (4.20) 
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This two-point function is scale- free and does not have the character of a canonical vacuum 
for a Rindler wedge. 

What conclusions can we draw from this? Previously, we argued heuristically that the 
SJ state in the Rindler wedge should be the Fulling-Rindler vacuum (to the extent that 
either is defined at all in the presence of the IR divergences). Now we have seen that a 
well-controlled limiting procedure gives a different result. As L gets large, the spacetime 
geometry approaches that of the Rindler wedge, as far as points that remain close to one 
corner of the diamond are concerned, but the SJ state approaches the ground state of a 
scalar field with refiecting boundary conditions at the corner. 

In fact, this "mirror behaviour" is already visible at the level of the SJ modes them- 
selves. The first set of modes fk (4.4) vanish on the two vertical lines at the spatial 
positions of the corners fk{x = ±V2L,t) = (recah that the corners are at x = 
in the coordinate system before translation), while the second set Qk also vanish on these 
vertical lines in the approximation /C — )• /Cq. The SJ conditions thus satisfy approximately 
the boundary conditions for two static mirrors, one at each corner of the diamond. How 
would such a two-mirror state appear near to the left corner? As the size of the diamond 
tended to infinity, one might expect the field in the left corner to become unaware of the 
right mirror, and this is consistent with our calculation above. On the other hand there 
remains the puzzle of where the left-hand mirror comes from in the limit. Its very existence 
selects a distinguished timelike direction, and since this direction can only be covariantly 
defined by reference to the right hand corner of the diamond, it seems difficult to avoid the 
conclusion that the presence of the right corner retains its influence no matter how large 
L becomes! 

It seems reasonable to attribute these counter-intuitive effects to our having set the 
mass to zero. As an aspect of its infrared pathology, the massless field might be able 
to sense the boundaries of the finite system, no matter how far away they are. If this 
explanation is correct, one would not expect to find the same mirror behavior for a massive 
scalar field, since the mass should shield it from such long range effects. It would also be 
interesting to study massless and massive fields in finite regions of Minkowski spacetime in 
3 + 1 dimensions. 

5 Comparison with the discrete SJ vacuum 

In this section, following [5, 6], we will apply the SJ formalism to the massless scalar field 
on a causal set that is well-approximated by the two-dimensional fiat causal diamond. In 
the case of non-zero mass, the second of these references has shown numerically that the 
mean of the discrete SJ two-point function approximates well the Wightman function of 
the continuum Minkowski vacuum. 

The discrete version of the SJ prescription can be interpreted in two ways. In the con- 
text of quantum gravity, causal sets are considered to be fundamental — more fundamental 
than continuum spacetimes which are just approximations to the true discrete physics [21]. 
On this view, the discrete SJ proposal is a starting point for building a theory of quantum 
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fields on the physical discrete substratum of spacetime and one can hope that it will give 
us clues about quantum gravity, just as quantum field theory in curved continuum space- 
time has done. From another viewpoint, the discrete formalism can be seen as simply a 
Lorentz-invariant discretisation of the continuum formalism and can therefore be used to 
test or extend the results of the continuum theory, in particular when analytic calculations 
are difficult. 

5.1 Causal sets and discrete propagators 

A causal set (or causet) is a set C together with an ordering relation ^ that satisfies the 
following conditions. It is reflexive: for sl\ X & C, X ^ X. It is antisymmetric: for all 
X,Y e C, X ^ X implies X = Y. It is transitive: for all X,Y, Z e C, X ^ Y ^ Z 
implies X ^ Z. And, it is locally finite: for all X,Y & C, \I{X, Y)\ < oo, where | • | denotes 
cardinality and I{X,Y) is the causal interval defined by I{X,Y) := {Z £ C\X :< Z < Y}. 
For more details on causal set theory the reader may refer to [21-23]. 

To produce a causal set Cm that is the discrete underpinning of (or approximation 
to, depending on the point of view) a given continuous spacetime M, we sprinkle points 
into M. A sprinkling generates a causal set from a given Lorentzian manifold by placing 
points at random in A4 via a Poisson process with "density" p. This produces a causal 
set whose elements are the sprinkled points and whose partial order relation is that of the 
manifold's causal relation restricted to the sprinkled points. The expected total number 
of elements in the causal set will be = pV_M. A causal set generated by sprinkling 
provides a discretisation of Ai which, unlike a regular lattice, is Lorentz invariant [24]. In 
the remainder of this paper, we will set p = 1 so that area will be measured directly by 
number of causet elements. Distance and area will thus be measured in natural units. 

The discrete SJ two-point function for a scalar field on a finite causal set [5, 6] is 
constructed using the same procedure as described above: from the retarded Green func- 
tion Gr, the Pauli-Jordan operator iA (a finite Hermitian matrix) is constructed and the 
Wightman function^ w is then obtained as the positive part of iA. In a finite causal set, 
this procedure is rigorously defined since everything is finite. It is of course necessary that 
the discrete analogue of the retarded Green function be known, which it is for massless [25] 
and massive [26] scalar fields on a sprinkling into two-dimensional Minkowski spacetime. 
The SJ two-point function is a complex-valued function of pairs of causal set elements, 
w : C xC — 7- C; equivalently it is a complex matrix w^^ := w{vi,Vj), where Vi^Vj G C. When 
C is obtained by sprinkling, each causal set element Vi corresponds to a point Xi in the em- 
bedding continuum spacetime. This allows us to directly compare the values of the discrete 
two-point function w^^ and those of the continuum Wightman function W^^ := W{Xi, Xj). 

There is numerical evidence [5] for a massive scalar field that at high sprinkling 
density, the SJ two-point function in the causet approximates the known Minkowski-space 
Wightman function. 

®We use a lower case w for the causet counterpart of the continuum function W. 
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Figure 4. An N = 2^^ = 2048 element sprinkling into a diamond Cl (with L = 2^ in natural 
units). The subregions corresponding to the centre («) and the corner (ii) are highlighted. 

We will extend this study to the massless case and compare with our results above for 
the continuum SJ vacuum. We will compare the mean of the massless discrete two-point 
function w with its continuum counterparts, using two separate methods. In the centre 
of the diamond, the continuum two-point function (4.16) is approximately a function of 
the geodesic distance only, in the limit of large L. We therefore plot the amplitudes w^^ 
against the proper time \d{Xi, Xj)\ and ask how well they agree with the continuum result. 
In the corner, the continuum W does not reduce, to a function of a single variable. In that 
case, we provide a "correlation plot" between the discrete two-point function and several 
continuum two-point functions (evaluated on the sprinkled points), so that the relative 
goodness of fit can be assessed. 

We restrict ourselves here to timelike related points; the analysis for spacelike related 
points is similar. Furthermore we only need to analyze the real parts of W and w. The 
imaginary parts add nothing new since they are given by the Pauli-Jordan function and 
tell us nothing about the quantum state. 

We work in a causal diamond M. = Cl and evaluate the discrete propagator on 
an N = 2^^ = 2048 element sprinkling into this diamond, which implies L = \/V /2 = 
= 25 in natural units. A typical sprinkling is shown in figure 4. Highlighted are the 
two subregions in which we shall sample the discrete S J two-point function. Each subregion 
occupies 8% of the area of the full diamond. 

5.2 The SJ vacuum in the centre 

The data for the centre is taken from a sample of 181 points in the square, among which 
there were 7599 timelike related pairs. As discussed in section 3 the two-point function of 
the massless scalar field in two-dimensional Minkowski spacetime is ill-defined owing to the 
infrared divergence, but there exists a one-dimensional family of "approximate Wightman 
functions" Wm,\ parameterized by an infrared scale or "cutoff" A. It is thus natural to 
compare our discrete function w with Wm,A; where we fix A = 0.02 from the relation (4.17) 



1.0 




0.8 



0.6 



0.4 



0.2 







2 



4 



6 



8 



10 



12 



Ml 



Figure 5. The real parts of the continuum two-point function Wm,\{X, X') (black line) with 
A — 0.02 and the discrete SJ two-point function Wij (blue scatter) in the centre of the finite 
diamond Cl with L — 2^ , plotted against the proper time \d\ for timelike separated events. 
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between A and L when L = 22. Then the real part of the continuum Wightman function 
we compare to is 



Figure 5 displays this function together with a scatter-plot of the discrete SJ amplitudes 



a deviation at larger values of proper time which, if real, can be attributed to 0{5/L) 
corrections, given that the centre region is still relatively large compared to L. From our 
previous analysis we know that Wm,\{X, X') approximates the continuum SJ state in the 
centre of the diamond, and so our data also supports the conclusion that the continuum 
and discrete SJ Wightman functions approximate each other. 

5.3 The SJ vacuum in the corner and in the full diamond 

For the corner, the type of plot we used for the centre is unsuitable because the continuum 
two-point functions we want to compare with depend on more variables than just the 
geodesic distance. Instead we plot the values of w directly against those of the continuum 
W with which we are comparing. More specifically, we use the coordinate values of the 
sprinkled points to calculate the values of a particular continuum function W{Xi, Xj) =: 
W^^ . We then plot a point on the graph whose vertical coordinate is W'^^ and whose 
horizontal coordinate is w^^ at the corresponding pair of elements of the causet. In this 
manner, we will assess the correlation between the data sets w^^ and W^^ for 4 different 
continuum two-point functions: the exact continuum SJ function Wsj^l (before Taylor 
expansion), the Minkowski function Wm,\ (3.6), the single mirror M^mirror (4.20), and the 
Rindler function Wr^\ (3.10). Both the Rindler and Minkowski W^- functions come with an 






Evidently, the fit is good, with a slight hint of 
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Figure 6. Correlation plots for the two-point functions in the corner of the causal diamond. The 
horizontal axis represents the causal set two-point function w'^ . The vertical axes represent, from 
left to right: the SJ, the Minkowski, the mirror, and the Rindler two-point functions. 

arbitrary parameter A, which shows up on the plots as an arbitrary vertical shift. We set 
this shift such that the intercept is zero. 

The corner subregion (figure 4) contains 181 points, which produce a sample of 11230 
pairs of timelike related points. The correlation plots for the real parts of the discrete and 
continuum propagators evaluated on this sample are shown in figure 6. Evidently, the plot 
of the exact SJ function exhibits a good fit with the causal set data, confirming again that 
the discrete and continuum formalisms agree. As expected, the correlation with the mirror 
two-point function is also high, implying that the ground state in the corner is indeed that 
of flat space in the presence of a mirror. (The slightly positive intercept in the causal set 
versus mirror plot, can plausibly be attributed to the e correction and to 0{S/L) effects, 
both of which would go away were the corner region made smaller.) As one would expect 
for the corner, the Minkowski and Rindler functions exhibit significantly worse correlations 
with the causet data-set. 

Turning to the full diamond, we use a smaller overall causal set with N = 256, which 
yields a sample of 16393 pairs of timelike related points. To the four comparison functions 
above we add a fifth: the reflecting box (mirrors at both corners) Wbox,L- Of these five con- 
tinuum two-point functions, one should expect that in addition to the SJ vacuum, only the 
reflecting box ground state would exhibit a reasonable degree of correlation with the causal 
set data. (As seen in equation (4.12) above, continuum SJ and reflecting-box are identical 
up to the error-term e(ti, v; u' , v'), which, however, can vary more appreciably now that we 
do not restrict ourselves to a small subregion of the diamond.) The correlation plots of 
flgure 7 conflrm this expectation, although the match with the discrete SJ function is not 
a sharp as in the case of the corner, perhaps reflecting the smaller overall sprinkling density. 

6 Concluding remarks 

We have determined the S J vacuum of a massless scalar field in the two-dimensional "causal 
diamond" spacetime and compared it with suitable reference states in the full Minkowski 
spacetime and in the Rindler wedge. Both of these reference "states", the "Minkowski 
vacuum" and the "Fulling-Rindler vacuum", suffer from infrared pathologies, but to the 
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Figure 7. Correlation plots for the two-point functions in the full causal diamond. The horizontal 
axis represents the causal set two-point function w^^ . The vertical axes represent, from left to right: 
the SJ, the Minkowski, the left mirror, the box (two mirrors), and the Rindler two-point functions. 



extent that they are meaningful we can compare them with the SJ state in the limit that 
the size of the diamond goes to infinity. 

When we look near the centre of the diamond, we find that the SJ Wightman function 
agrees with the "Minkowski vacuum" , just as one might have expected. However, when we 
look in the corner of the diamond we do not see something having the form of a "Fulling- 
Rindler vacuum". Instead we recognize the flat-space vacuum in the presence of a static 
mirror at that corner. This is confirmed by the numerical results. Thus, the continuum 
calculations in section 2 giving the SJ state in the Rindler wedge do not agree with the 
limiting procedure of constructing the SJ state in the finite diamond and letting the size 
of the diamond tend to infinity, whilst keeping the arguments of the two-point function at 
fixed locations with respect to the corner. It is important to understand the reason for 
this ambiguity which threatens the proposal of the SJ state as the distinguished ground 
state of a region: such a ground state, if it is defined at all, should be unique and be 
independent of any (physically sensible) limiting procedure. Of course, one cannot really 
speak of a disagreement between two functions, one of which is ill-defined, but the fact 
remains that our limiting procedure yields very different results for the centre versus the 
corner of the diamond. It seems likely that these disagreements and ambiguities stem from 
the infrared divergences of the two-dimensional massless theory, and if one were to work 
instead with a massive field, ^ the SJ state for the wedge would be unique and in agreement 
with both the Fulling-Rindler vacuum and the limit of the diamond's SJ state. Some of 
these questions will be investigated in future work, but preliminary numerical results for 
the massive scalar field in the sprinkled diamond indicate that the discrete SJ function 
does fit the Fulling-Rindler vacuum better than the state with a mirror present. 

If the foregoing expectations are born out, then one conclusion would be that the 
SJ state for the massive field is singular on the boundary of the diamond, which is the 
behaviour one would expect for a pure state in a bounded region. (The S J state of a region 
is pure by construction.) Indeed, the highly entangled nature of quantum states possessing 
the local structure of the Minkowski vacuum means that their restrictions to any spatially 
bounded portion of Minkowski spacetime will be highly mixed and far from pure. 

We end by considering again the "tug of war" that Fulling identifies in the theory of 
quantum fields. The split of spacetime into space plus time that Fulling assumes to be 

''One might also consider working with the gradient of the field. 
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inherent in the quantum aspect of the field theory is actually an artefact of the choice to 
conceive of a quantum field as if it were the canonical quantisation of a classical Hamil- 
tonian system. In a canonical approach, defining the Hamiltonian tends to demand the 
foliation of spacetime into spacelike hypersurfaces. However, there is an alternative, long 
ago identified by Dirac as more fundamental because it is essentially relativistic: the quan- 
tum analogue of the classical Lagrangian approach, namely the path integral [27]. Like the 
path integral, the SJ vacuum is also relativistic in character, depending, as it does, on the 
geometry of the whole spacetime region via the causal structure encoded in the retarded 
Green function. Its construction is inherently covariant, and at no point is there a need 
to introduce mode-functions defined on a hypersurface, except as dictated by calculational 
convenience. The formulation of the quantum theory along these lines bears no resem- 
blance to the "canonical quantisation" process: one does not solve the classical equations 
of motion or identify canonically conjugate variables or promote them to operators satis- 
fying canonical commutation relations. Such a formulation seems much more compatible 
with a path integral approach to quantum theory, and indeed the SJ proposal serves as the 
starting point for the construction of a histories-based formulation of quantum field theory 
on a causal set [3] which admits a natural generalization to interacting scalar fields and 
takes us one step further towards a quantum theory of causal sets. 
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A Corrections to the SJ two-point function 

In evaluating the second sum (4.9) of the continuum SJ two-point function (4.6), we made 
the approximation /C — )• /Co- For a given pair of spacetime points, this will induce an error 
in the two-point function, which (see (4.10) and (4.6)) is given by 



e{u,v;u,v') = ^ 



n=l 



:9k„ (u, v)gl^ (n , v') - ^ na gfco.n (^> v)gl^ ^ {u , v) 



k \\m 2^'-"\-' -^-^Kn^" ' " ^ n, 
'^n \ \ykn 1 1 '^0,11 I \ykQ , 

(A.l) 

where kn and ko,n denote the n*'' terms in /C and /Co, respectively. The contributions to the 
right hand side come mostly from long wavelength (small n) terms, since the approximation 
/C — >• /Co becomes increasingly accurate for large n (see figure 3). This means that we 
should expect e to be constant over small subregions of the diamond. We will first test 
this expectation numerically, restricting ourselves for simplicity to timelike related pairs of 
points 

To estimate the mean and the variation e{u,v;u' ,v') over different pairs of spacetime 
points in a subregion associated with the centre (i) or corner (ii) of the diamond, we 
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Figure 8. The mean and standard deviation of e{u,v;u' ,v') in the centre (left) and in the corner 
(right). The vertical axis corresponds to e{u,v;u' ,v') and the horizontal axis is — logio(Kiub/l^)- 

evaluated e{u,v\u' ,v') on a random sample P of pairs of timelike related points within 
that region. We evaluated (A.l) by truncating the sums on the right hand side at a stage 
large enough for the sum to have converged sufficiently. We then calculated the mean 
value of e{u,v;u',v') and its standard deviation on the sample P. We present the results 
for the real part of e{u,v;u',v') below, since we are mainly interested in the real part of 
the Wightman function W. (The imaginary part of W is proportional to A and therefore 
known exactly. Moreover, the imaginary part of e{u,v;u' ,v') was consistent with zero in 
all the regions we investigated numerically.) 

Look at subregions in the centre and corner of the form depicted in figure 4: a square 
in the centre, a triangle in the corner. Fix the linear dimension of the subregion D L 
and denote its spacetime volume by V^ub- Increase the size L of the full diamond while 
keeping D fixed, thereby decreasing the volume ratio V^ub/^- The mean and standard 
deviation of e{u,v;u',v') obtained in this way are shown in figure 8 for different values of 
Vsuh/y ranging from 0(10"^) to 0(10"^). These results were obtained by truncating the 
sum (A.l) at n = 50, which provides sufficient accuracy. We observe that the standard 
deviation in e{u,v;u' ,v') indeed quickly becomes negligible as Kjub/^ is decreased. The 
mean of e{u,v;u' ,v') tends to a constant value in the centre given by ecentre = —0.0627 
and it vanishes in the corneri ^corner — 0. Notice that these results are unchanged under 
a simultaneous rescaling of L and D: the mean and standard deviation of €{u,v;u',v') 
depend only on the ratio Vsuh/y- 

It is worth noting that the asymptotic values for the mean of e seen at large L above 
agree with the values of e{u, v; u', v') that we obtain if we simply evaluate the infinite sum 
on a pair of coincident points in the centre of the diamond, {u,v) = {u',v') = (0,0), or in 
the corner {u,v) = {u',v') = {—L,L). In the centre, the sum (A.l) reduces to 




(A.2) 
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where x„, is the n positive solution to tan(x) = 2x. This sum can be evaluated to 
arbitrary precision using numerical solutions for x„, and it tends to e(0, 0; 0, 0) = —0.0627, 
corresponding to the horizontal asymptote in the centre plot of figure 8. In the corner, 
both terms in (A.l) vanish because the modes (4.4) are identically zero at {—L, L) for all 
kn G JC and A;o,n G A^o- It follows that e(— L, L; —L, L) = 0, corresponding to the horizontal 
asymptote in the corner plot of figure 8. 
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